Observation of Parity Violation in the Q — > AK Decay 
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Abstract 

The a decay parameter in the process Vt~~ — > AK~ has been measured from a sample 
of 4.50 million unpolarized f2~ decays recorded by the HyperCP (E871) experiment 
at Fermilab and found to be [1.78±0.19(stat)±0.16(syst)] xlCT 2 . This is the first 
unambiguous evidence for a nonzero a decay parameter, and hence parity violation, 
in the Q~ — > AK~ decay. 
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Our knowledge of the Q~ hyperon and its decays remains incomplete, despite 
its long and illustrious role in particle physics. Its spin and parity have not been 
firmly established, 1 and it alone among the hyperons has yet to exhibit parity 
violation in its two-body weak decays. The Particle Data Group (PDG) values 
of the a decay parameters of the three such decays, fi~ — > AK~, fi~ — > S _ 7r°, 
and n~ -> S°7T-, respectively -0.026±0.023, +0.05±0.21, and +0.09±0.14 
[2], are consistent with zero, where a is a measure of the interference between 
the P- and D-w&ve final-state amplitudes: 



2Re(P*D) 

a= \pJTW () 



A nonzero value of a is manifest evidence of parity violation. In contrast, all 
other hyperons have been shown to have nonzero a decay parameters. The 
smallest are those of the S + — > im + and the Tr — > rm~ decays, both of which 
are 0.068; the largest is almost unity: a = —0.980 in S + — > pir decays [2]. The 
two-body nonleptonic Q~ decays are expected to be nearly parity conserving 
[3], and hence predominantly P wave, implying a small a decay parameter, 
which is consistent with the experimental results. 

Recently, we have reported evidence of parity violation in an analysis of 0.96 
million Q~ — > AK~ decays taken in the 1997 Fermilab fixed-target running 
period, yielding a n = [2.07±0.51(stat)±0.81(syst)] x 10~ 2 [4]. (Throughout 
this Letter an will refer only to the AK~ decay mode of the We report 
here another measurement of «o using 4.50 million events taken during the 
1999 Fermilab fixed-target running period. 

The experiment was mounted in the Meson Center beam line at Fermilab us- 
ing an apparatus [5] built to search for CP violation in hyperon decays (see 
Fig. 1). A negatively charged secondary beam with an average momentum 
of 160GeV/c was produced by steering an 800GeV/c proton beam onto a 
60mm long, 2x2 mm 2 wide, Cu target. The target was followed by a curved 
collimator embedded in a dipole magnet ("Hyperon Magnet"). The f2~'s were 
produced at an average angle of 0°, ensuring that their polarization was zero. 
The secondary beam exited the collimator upward at 19.51 mrad relative to 
the incident proton beam direction. A 13m long evacuated pipe ("Vacuum 
Decay Region") immediately followed the collimator exit. After the Vacuum 
Decay Region was a high-rate magnetic spectrometer employing nine multi- 
wire proportional chambers (MWPCs), four in front of a pair of dipole magnets 
("Analyzing Magnets"), and five behind. Particles with the same sign as the 



* Corresponding author. Tel.: 1-434-982-5364, fax: 1-434-982-5375; E-mail address: 
craigdukes@virginia.edu (E.C. Dukes). 

1 The Q~ spin has not yet been determined, but measurements have ruled out J = \ 
and are consistent with the quark- model prediction of J = |; see [1]. Throughout 
this Letter we assume that the is spin-|. 
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Fig. 1. Plan view of the HyperCP spectrometer. 

secondary beam were deflected by the Analyzing Magnets to the beam-left 
side of the apparatus, and those with the opposite sign to beam-right. The 
highly redundant tracking system facilitated very high track-reconstruction 
efficiencies. The trigger required the coincidence of at least one hit counter in 
each of the Same-Sign (SS) and Opposite-Sign (OS) hodoscopes situated on 
either side of the secondary beam (the LR subtrigger), along with an energy 
deposit of at least ^40 GeV in the hadronic calorimeter. This energy threshold 
was well below the 60 GeV of the lowest-energy protons from f2~ decays, all 
of which entered the calorimeter, and above the energy where the calorimeter 
efficiency plateaued at ~99%. Since there was a high probability that both the 
K~ and the n~ would hit the SS hodoscope and since the OS hodoscope had 
two layers of counters, the efficiency of the LR subtrigger was extremely high 
(^99.5%). Events that satisfied the trigger were written to magnetic tape by 
a high-rate data acquisition system [6]. 

The analysis reported here is from data taken with the negative-polarity sec- 
ondary beam. The 29 billion recorded events were initially reconstructed and 
separated according to event type using loose event-selection cuts. This left 
a total of 56 million candidate events. The raw event information was pre- 
served at this (as well as every subsequent) stage. Final event-selection crite- 
ria were applied after careful study and were tuned to maximize the signal- 
to-background ratio. The most important requirements were that: (1) the 
X 2 /df for a geometric fit to the decay topology be less than 2.5; (2) the 
distance-of-closest-approach for the tracks forming the A and Q~ decay ver- 
tices be less than 4 mm; (3) the x and y separations from the target cen- 
ter of the extrapolated Q~ trajectory satisfy the inequality (i/2.0mm) 2 + 
(«//2.2mm) 2 < 1.0; (4) both the f2~ and the A decay vertices lie at least 
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0.28 m (0.32 m) downstream (upstream) of the entrance (exit) of the Vacuum 
Decay Region, and that the Q~ vertex precede that of the A; (5) the pir~7T~ 
(7r + 7r~7r~) invariant mass be greater than 1.355 GeV/c 2 (0.520 GeV/c 2 ), in or- 
der to eliminate 5~ — > An~ — > pn^n^ (K~ — > 7r + 7r~7r~) decays; (6) the pn~ 
and pir~K~ invariant masses be respectively within ±4.0 MeV/c 2 (4.3 a) and 
±8.0MeV/c 2 (5.0 a) of the true A and Q~ masses; and (7) no particle have 
momentum less than 12GeV/c. After all these cuts the number of events re- 
maining was 4.735 million. Monte Carlo simulation indicated that 55.3% of 
Q~ — > AK~ — > pn~K~ decays for which the Vt~ exited the collimator passed 
these cuts. The cuts that eliminated the greatest numbers of signal events 
were the pn^n^ invariant mass and the Q~ vertex requirements. 

Figure 2 shows the pn~K~ and p7i~ invariant-mass distributions after event 
selection cuts. The background-to-signal ratio, determined using a double- 
Gaussian plus second-degree polynomial fit to the invariant-mass distribution, 
is (0.33±0.03)% in the region within ±5.0 a of the Q~ mass. The background 
under the p7i~ mass peak is less than half this. Dominant backgrounds were 
misreconstructed H~ — > An~ — > pn^7i~ decays and Q~ — > H°7r _ — > A7r°7r~ — > 
p7r~7r~77 decays. 
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Fig. 2. The pir~ K~ (left) andp7r~ (right) invariant-mass distributions, after all cuts 
except the respective mass cuts. Arrows delimit the accepted mass regions. 



The Vt~ alpha parameter was measured through the asymmetry in the A — > p-n~ 
decay distribution. In the decay of an unpolarized fl~ to a A and a K~, the 
A is produced in a helicity state with its helicity given by [7]. Hence the 
decay distribution of the proton in that A rest frame in which the A direction 
in the Q~ rest frame defines the polar axis — the Lambda Helicity Frame 
shown in Fig. 3 — is given by 

— (1 + a n a A cos9), (2) 



d cos 9 2 
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where 9 is the polar angle of the proton and «a is the alpha decay parameter 
in A — > pir~ . Since the A decay direction in the Q~ rest frame changes from 




Fig. 3. The Lambda Helicity Frame. 

event to event, so too does the polar axis of the Lambda Helicity Frame along 
which the A polarization must lie: knowledge of the direction of the putative 
A polarization is of enormous importance as it greatly minimizes biases. The 
analysis "locks in" to the changing direction of the A polarization. Biases, on 
the other hand, such as uncorrected detector inefficiencies, are fixed in the 
laboratory frame. Hence the Lambda Helicity Frame analysis acts much like 
a lock-in amplifier, except that it locks into a known direction rather than a 
known frequency. 

The proton cos 6? acceptance was measured and corrected for using a hybrid 
Monte Carlo (HMC) technique that has been used in many similar such mea- 
surements [8]. Monte Carlo events were generated by taking all parameters 
from real events except the proton and pion direction in the rest frame of 
the A. An isotropic A — > pir~ decay was generated, and the proton and pion 
were boosted back into the laboratory frame using the real A momentum. 
Their trajectories were then traced through the apparatus, with the detector 
responses simulated where appropriate (using measured efficiencies), and all 
MWPC wire hits not associated with the real proton and pion tracks kept. The 
simulation included multiple scattering and slope-dependent multiple-wire hit 
probabilities which were tuned to match data. The HMC simulated the data 
extremely well, as is evident by the small x 2 /df in the matching of the real and 
HMC proton distributions in the Lambda Helicity Frame (see discussion be- 
low). Real and HMC distributions of proton and pion tracks at various places 
along the spectrometer also showed excellent agreement. The HMC proton 
and pion tracks, in conjunction with the real kaon, were required to satisfy 
the trigger requirements, and were reconstructed by the standard track-finding 
program, with the same cuts applied to all parameters formed from them that 
were applied to the real events. Ten accepted HMC events were used for each 
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real event. If over 300 generated HMC events were required to get the ten, 
then both the real and associated HMC events were discarded; this eliminated 
regions of low acceptance and reduced the computer time needed for the anal- 
ysis. It eliminated 4.9% of the events. Increasing the upper limit beyond 300 
events had no effect on the result. 

Since the HMC events were generated with a uniform proton cos 9 distribution, 
each accepted HMC event was then weighted by 

W= 1 + Scos0 f (3 ) 

where S is the slope (to be determined) of the proton cos^ distribution and 6f 
and 6 r are, respectively, the HMC ("fake") and real proton polar angles in the 
Lambda Helicity Frame. Note that in the absence of a background correction 
S = a^a\. The numerator in Eq. (3) in effect polarizes the HMC sample, while 
the denominator removes the polarization bias accrued from using parameters 
from real polarized A decays. To facilitate handling the unknown slope S, 
the weights, binned in cos#j, were approximated by the polynomial series 
expansion 

W « (l + Scos9 f )[l - Scos6 r + (Scos6 r ) 2 + (S cos6 r ) 10 }. (4) 

The polynomial coefficients, which depend only on cos#j and cos# r , were 
summed, and then S was extracted by minimizing the x 2 difference between 
the real and weighted HMC proton cos 9 distributions. The error was deter- 
mined by finding the variation in S needed to increase x 2 by one. It includes 
the uncertainty in the acceptance as determined by the HMC events. 

The analysis procedure was extensively checked by Monte Carlo simulation. 
Monte Carlo fi~ — > AK~ — > pn~K~ events were simulated using the real 
hodoscope, MWPC, and calorimeter efficiencies, and required to pass the same 
cuts as the real data. These were analyzed by the HMC analysis code. The 
input and extracted values of a^ctA were found to be consistent over a wide 
range of «o input values, with an average difference of (0.017±0.042) x 10~ 2 . 
As a cross-check, 78 000 S~ — > An~ — > pir~iT~ decays available from the same 
dataset were analyzed using exactly the same analysis program, with selection 
criteria tuned for the S~ decay. The fit to the proton cos 6 distribution was 
good, with x 2 /df = 14/19. The correct sign of a~,ct\ was found, which is 
opposite the sign of our value of a^a\, and the magnitudes of the measured 
and PDG values of a^a\ were consistent within the statistical errors. 

A total of 4 504 896 real events were analyzed by the method described above. 
The real and weighted HMC proton cos# distributions are shown in Fig. 4, 
and the differences between the real and HMC proton cos6 distributions, 
weighted and unweighted, are shown in Fig. 5. The nonisotropic nature of the 
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Fig. 4. The real (lines) and weighted HMC (points) proton cos 9 distributions. The 
total number of HMC events has been scaled down by a factor of 10. 

real proton cos 6 distribution, compared to the isotropically generated HMC 
distribution, is clear from the top plot of Fig. 5. It is unambiguous evidence of 
a nonzero a decay parameter. The bottom plot shows the same comparison, 
except that the HMC events have been weighted by the best-fit value of S. 
The extracted slope of the proton cos# distribution is S = (1.16±0.12)xl0~ 2 
with x 2 / df = 23/19, where the error is statistical. 
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Fig. 5. The relative differences between the real (N r ) and HMC (Nf ) proton cos 9 
distributions for unweighted (top) and weighted (bottom) HMC events. The total 
number of HMC events has been scaled down by a factor of 10. 
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To extract asj«A from the proton cos# slope, the small background contribu- 
tion to S was subtracted. To estimate the proton cos 9 slope from the back- 
ground events the same analysis procedure was performed on five sideband 
regions, three below and two above the Q~ mass region. The average side- 
band proton cos 8 slope was found to be S^b = (7.2±3.0)xl0 -2 , with average 
X 2 /df = 19/19. No mass dependence of was apparent. The contribution 
to S of the background under the mass peak was corrected for by subtracting 
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the appropriate fraction of S s b from S, giving a^a\ = [1.14±0.12(stat)]xl0 
Note that this correction is only a 1.7% effect. 
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Fig. 6. The value of aaa\ vs. the momentum. 



The stability of the result was studied as a function of several parameters. 
The value of a^a\ was independent of the Q~ momentum, as shown in Fig. 6, 
and there was no dependence on the z location of the Q~ decay vertex. The 
non-background-subtracted slope S, measured on a run-by-run basis for all 
175 runs in the dataset, shows no evidence of a temporal dependence (see 
Fig. 7). 
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Fig. 7. Run-by-run values of the proton slope S. The dashed line represents the 
best-fit value. 




Systematic errors were small because of the high efficiencies of the spectrome- 
ter elements and because of the power of the Lambda Helicity Frame analysis. 
The dominant systematic errors are listed in Table 1. The effects of uncertain- 
ties in detector inefficiencies — MWPCs, trigger hodoscopes, and hadronic 
calorimeter — on «[)«a were found to be negligible. No statistically signifi- 
cant difference in S was found between using perfect and measured detector 
efficiencies when simulating the HMC proton and pion. The combined effect of 
the uncertainties in the fields of the Analyzing Magnets, 5.5 G, was also negligi- 
ble. A small fraction of the daughter 7r~'s and K~ J s decayed before exiting the 
apparatus. (Approximately 0.7% of the 7r~'s decayed before the last MWPC.) 
The effect of such decays on a^a\ was studied using Monte Carlo events and 
data and found to be negligible. The error in the background subtraction was 
estimated by assuming that the error in the average sideband slope S s ^ was 
equal to the average sideband slope, 5S s b = 7.2xl0 -2 , and using a 25% error 
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in the background-to-signal ratio, both very conservative assumptions. It too 
is negligible. 

Table 1. Systematic errors. 



Source Error (10 -2 ) 

Event selection cut variations 0.088 

Validation of analysis code 0.042 

Background subtraction uncertainty 0.024 

Detector inefficiency uncertainties O.OfO 

Analyzing Magnets field uncertainties 0.006 



The largest systematic uncertainty was the sensitivity of the measurement 
to the values of the cuts used to define the data sample. The most impor- 
tant were the cuts on the piT~ and piT~K~ invariant masses and, less im- 
portantly, the 12GeV/c minimum momentum cut. The effect of changes in 
these cut values was 0.088xlCT 2 . The total systematic error, including the 
upper limit in the uncertainty of the MC validation of the analysis program 
(0.042x 10 -2 ), is O.lOx 1CT 2 . This is a factor of five reduction in the systematic 
error of 0.52 xlO -2 reported in the analysis of the 1997 data [4]; most of the im- 
provement comes from incorporating the measured detector and track-finding 
inefficiencies into the HMC simulation in this analysis. 

To summarize, we find from a sample of 4504 896 Q~ — > AK~ — > pn~ K~ 
decays the value ooOa = [1.14±0.12(stat)±0.10(syst)] x 10~ 2 . Using oa = 
0.642±0.013 [2], a n is found to be [1.78±0.19(stat)±0.16(syst)] x lO" 2 , where 
the contribution of the uncertainty in the value of a\ to the systematic error 
is negligible. Our measurement represents a factor of nine improvement in 
precision over the current PDG value. It is 1.9 a from the PDG average of 
(— 2.6±2.3) x 10~ 2 , and opposite in sign. This measurement is consistent with 
the recent result we reported [4] from an independent analysis of data taken in 
the 1997 fixed-target running period, but with a factor of four smaller error. 
With a magnitude that is 7.2 a from zero, it represents unambiguous evidence 
of parity violation in the Q~ — > AK~ decay. As predicted, «jj is small; indeed 
it is the smallest of all the a parameters that have been measured in the 
two-body weak decays of hyperons. 
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